en 



43 

Of 



Is the Weibull distribution really suited for wind 
statistics modeling and wind power evaluation? 

P. Drobinski 1 and C. Coulais 2 

1 Laboratoire de Meteorologie Dynamique - Institut Pierre Simon Laplace, CNRS 
and Ecole Polytechnique, Palaiseau, France 

2 CEA Saclay/DSM/IRAMIS/SPEC/SPHYNX, 91191 CEDEX Gif-sur-Yvette, 
France 

E-mail: philippe . drobinski@lmd . polytechnique . f r 



Abstract. Wind speed statistics is generally modeled using the Weibull distribution. 
This distribution is convenient since it fully characterizes analytically with only two 
parameters (the shape and scale parameters) the shape of distribution and the different 
moments of the wind speed (mean, standard deviation, skewness and kurtosis). This 
^ . distribution is broadly used in the wind energy sector to produce maps of wind energy 

^->| potential. However, the Weibull distribution is based on empirical rather than physical 

justification and might display strong limitations for its applications. The philosophy 
of this article is based on the modeling of the wind components instead of the wind 
speed itself. This provides more physical insights on the validity domain of the Weibull 
distribution as a possible relevant model for wind statistics and the quantification of 
the error made by using such a distribution. We thereby propose alternative expres- 
sions of more suited wind speed distribution. 

Keywords: Surface wind, Wind statistics, Wind energy, Weibull distribution, Rayleigh 
distribution, Rice distribution, Supcrstatistics. 
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1. Introduction 



Understanding and modeling wind speed statistics is a key for a better understanding 
of atmospheric turbulence and diffusion and for use in such practical applications as 
air quality and pollution transport modeling, estimation of wind loads on buildings, 
prediction of atmospheric or space probe and missile trajectory and wind power analysis. 
In the context of wind energy production, it is necessary to characterize the wind in a 
planned wind turbine farm location, in order to calculate the optimal cut-in and cut-out 
speed and the likely power output of a given wind turbine (Drobinski 2012). This is 
done using the distribution of wind speed M. In particular, a key parameter to describe 
wind power is the average wind power produced by a wind turbine, 

P w = ^pAC p M 3 K ej (1) 
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where A is the area of the flow, p is the density of the fluid, C p is the coefficient of power 
of the wind turbine which measures how efficiently the wind turbine converts the energy 
in the wind into electricity, K e is the energy pattern, and M 3 is the third moment of 
the distribution. In order to describe the wind energy potential independently of the 
turbine shapes, observed histograms are fitted to the Weibull distribution, defined as 
follows, 



where M is the wind speed, k > is the shape parameter and c > is the scale parameter 
of the distribution. In this case, M 3 = c 3 T (1 + 3/k) and K e = T (1 + 3/Jfc) /T 3 (1 + l/k), 
where T is the Gamma function. The Weibull distribution is thus an extremely used 
paradigm to model wind statistics (e.g. Justus et al. 1976, 1978, Seguro and Lambert 
2000, Cook 2001, Weisser 2003, Ramirez and Carta 2005) because it can be easily fit 
after estimating the parameters k and c and it produces analytical expressions for the 
different moments and the wind power. However, its use is purely empirical and there 
is a lack of physical background justifying the use of the Weibull distribution to model 
wind statistics. Besides, the validity of the description of wind statistics by a Weibull 
distribution is not clearly established. Therefore, the estimated wind potential might 
be systematically mis-estimated. 

In this paper, we use meteorological data to show that, in some locations with 
particularly high wind potential but wind anisotropy, wind statistics is not well described 
by the Weibull distribution, thus leading to a systematic under-estimation of wind 
power, up to 12%. We then provide physical insights for the origin of such a discrepancy. 
To this end, we consider wind velocity as a bivariate distribution of the two orthogonal 
(e.g. north-south and east- west) wind components and we compute it analytically 
in the simple case of normal distributions. This allows us to derive a Rayleigh-Rice 
super-statistical distribution, which is better suited than the Weibull distribution to 
describe real data in locations where the wind is not isotropic. Finally, we generalize 
this approach to the case of isotropic wind with non-Gaussian wind components, where 
the Weibull distribution better describes the meteorological data. 

2. Observations 

In this study, we use wind measurements collected between December, 28 s * 2006 and 
March 20 th 2008 at 10 m height from 4 surface weather stations at different locations 
in France (Fig. d|). Among the 4 stations, 3 are operational weather stations from 
the French national weather service Meteo-France and one is the research observatory 
(Haeffelin et al. 2005). While the wind components from the 3 operational weather 
stations are recorded every three hours, the SIRTA observatory data is recorded every 
minute. The SIRTA observatory is located at Palaiseau, 20 km south of Paris, France. 
The site is in a suburban area. We use the wind measurements on a tower surrounded 
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Figure 1. Region of interest with the shaded area indicating a topography higher than 
500 m. Locations of the surface weather stations of the SIRTA research observatory, 
25 km South of Paris, and from the Meteo-France network at Lyon, Montelimar and 
Orange are indicated with a " x " . 



by different types of surface, i.e. close and distant forest, buildings and an open field 
sector (Drobinski et al. 2006, Barthlott et al. 2007, Fesquet et al. 2009). 

The Meteo-France stations are located in a much more complex environment with 
elevated mountain ridges. In this region, frequent channelled wind can persist for several 
days. The strongest and most frequent channeled valley wind is the mistral, which blows 
from the north/north- west in the Rhone valley, separating the French Alps (highest 
elevation 4807 m) to the east from the Massif Central (highest elevation 1885 m) to the 
west by a gap 200 km long and 60 km wide (Fig. [T|) (Drobinski et al. 2005). It occurs 
when a synoptic northerly flow impinges on the Alpine range. As the flow experiences 
channeling, it is substantially accelerated and can extend offshore over horizontal ranges 
exceeding few hundreds of kilometers (Salameh et al. 2007, Lebeaupin-Brossier and 
Drobinski 2009). The mistral occurs all year long but exhibits a seasonal variability 
either in speed and direction, or in its spatial distribution (Mayencon 1982, Orieux and 
Pouget 1984, Guenard et al. 2005, Guenard et al. 2006). During summer, the mistral 
shares its occurrence with a northerly land breeze and southerly sea-breeze (e.g. Bastin 
and Drobinski, 2005, 2006; Drobinski et al., 2006, 2007), which can also be channelled 
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Figure 2. Probability density functions (PDF) of wind speed (blue) in m s _1 and a 
fit by a Weibull function (dashed black) from observations made at SIRTA (a), Lyon 
(a), Montelimar (c) and Orange (d). 



in the nearby valleys (Bastin et al., 2005) or interact with the mistral (Bastin et al., 
2006). In such region, accounting for such persistent wind systems for modeling the 
wind speed statistics is thus mandatory. 

Figure [2] shows the probability density functions (PDF) of the wind speed at 
the SIRTA observatory, Lyon, Montelimar and Orange and their fits by a Weibull 
distribution. First, while the Weibull distribution is an accurate statistical model for 
the wind speed distribution at SIRTA (Fig. EK), it poorly describes wind distributions 
measured at the 3 weather stations located in complex terrain (Fig. |2b-d). Table Q] 
compares the wind power computed from the wind speed observations and that is 
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Table 1. Comparison of the power computed from the measured data and out 
of the different fits (Weibull, Elliptical and Raylcigh-Ricc distributions; see text for 
definitions) using typical values Cp = 0.45, p = 1.2 kg/m 3 and A = 7.9 x 10 3 m 2 . The 
quantities k and c are the estimated Wcilbull distribution parameters, a u and a v the 
estimated standard deviations of the 2 wind components of the Elliptical distribution 
and a is the estimated regime parameter (see text for definition). 



Stations 




Powei 


• x(10 2 kW) 






Data 


Weibull 


Elliptic 


Raylcigh-Ricc 






[k;c] 




[a] 


SIRTA 


1.63 


1.59 


1.48 


1.13 






[1.71;3.58] 


[2.64;2.64] 


[0.05] 


Lyon 


2.26 


2.22 


2.17 


2.26 






[1.76;4.11] 


[2.33;3.54] 


[0.297] 


Montclimar 


2.75 


2.55 


2.53 


2.74 






[1.82;4.48] 


[2.65;3.73] 


[0.230] 


Orange 


7.71 


6.78 


6.75 


7.59 






[1.45;5.69] 


[2.93;5.53] 


[0.510] 



estimated from the Weibull distribution best fit using Eq. (CQ). The parameters k and c 
are also indicated. At the SIRTA observatory, the wind power estimated from Eq. (pQ) 
is accurate up to 2.5 %. It is obtained for a value k = 1.71. The Weibull distribution 
describes fairly well the wind distribution at SIRTA, which we stress, is located in a fairly 
flat terrain. Yet, the difference between the data and the fit is not random fluctuations: 
there are small systematic differences occurring at large (M > 8 m.s -1 ) and at low 
(M < 2 m.s -1 ) wind speed values. 

In contrast, at Montelimar and Orange, Weibull fit quality drops (Fig. [2b-d), and 
systematically underestimates the right tail of the distributions. As a result, the power 
is systematically underestimated by 7.5 % and 12.8 % (Tabled]). To the best of our 
knowledge, this systematic discrepancy has not been reported in the existing literature, 
although it is a major issue for wind power resource evaluation. At Lyon, Montelimar 
and Orange, the observed distributions display complex shapes instead of a smooth 
and regular shape as would be expected when fitting with a Weibull distribution. As 
opposed to the regular shape of wind speed and the surrounding topography at SIRTA, 
the complexity of the distributions in the Rhone valley might originate in the complexity 
of the surrounding topography. For a deeper insight of the differences between observed 
wind speed distributions and the commonly used Weibull distribution, the distribution of 
the wind components are now investigated (Fig. [3]). At the SIRTA observatory (Fig.[3h/), 
the zonal and meridional wind components have the same symmetrical distributions with 
zero mean and equal width: the wind is isotropic. In addition, the distributions of wind 
components have larger tails than Gaussian. 

At Lyon (Fig. [3b), at the entrance of the Rhone valley, the PDF of the zonal wind 
component is still symmetrical with almost zero mean and has a Gaussian shape. The 
PDF of the meridional wind component is fairly symmetrical for values below 5 m s _1 but 
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Figure 3. Distributions of north-south (red) and East- West (blue) wind components 
in m s _1 and a Gaussian function (dashed black) from observations made at (a) SIRTA, 
(b) Lyon, (c) Montelimar and (d) Orange. Insets : Joint distributions of wind 
components. 

is skewed for the strongest values. Also the width of the meridional wind PDF is larger 
than that of the zonal wind. At Montelimar and Orange, none of the wind component 
PDFs look like Gaussian distribution and the meridional wind PDF become significantly 
larger than the PDF of the zonal wind. This shows evidence of the channelled nature 
of the wind by the north-south oriented Rhone valley (Salameh et al. 2009). For low 
winds, the PDFs of the zonal and meridional components are more similar, evidencing 
a less efficient impact of the valley flanks to orientate the wind. 

The two wind components may however be correlated. The insets of Fig. [3] show the 
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joint distribution of the two wind components at SIRTA observatory, Lyon, Montelimar 
and Orange. At the SIRTA observatory, the rather flat terrain produces an isotropic 
wind statistics (i.e. the wind blows in all directions) with no correlation between the 
two components (i.e. the joint distribution displays a circular shape). Conversely, the 
joint distributions at Lyon, Montelimar and Orange display strongly correlated wind 
components in specific directions. The main directions correspond to the north-south 
orientation of the Rhone valley. 

To conclude on the observations, it appears that at SIRTA and Lyon, while the 
Weibull function fits better wind speed distributions, the wind components are less 
anisotropic. Therefore, it seems that isotropy is crucial for the Weibull description to 
be valid. In the following, despite the fact that none of the observed distributions is 
Gaussian, we shall focus on anisotropic bivariate distribution to understand the effect 
of anisotropy on the resulting wind speed distribution. 

3. From Gaussian to super-statistical Gaussian wind components 

3.1. The Rayleigh distribution 

The very first approach to model anisotropy is to consider a bivariate distribution 
of Gaussian wind components with means \x u and [i v , and variances a u and a v and 
finite cross-correlation p. After a change of variable, u — > u', v — > v', the resulting 
joint distribution is a bivariate elliptical distribution with variances a' u and a' v without 
correlation. In the case p! u = fi' v = (Crutcher and Baer 1962), 

Applying the usual transformation from Cartesian to polar coordinates (M,0), and 
integrating over the angle 0, the joint PDF for M and <fi is (Chew and Boyce 1962): 

p M (M; o», O = ex P H M2 ) J o (bM 2 ) , (4) 

with a = (c' 2 +a' 2 )/(2a' u a' v ) 2 and b = (al—a 2 )/(2a u a v ) 2 , and where I Q (x) is the modified 
Bessel function of the first kind and zero order. When the variance are equal 6 = 0, 
so I (bM 2 ) = 1 and Eq. (jl]) is the Weibull distribution with parameter k = 2, which 
is the Rayleigh distribution. When comparing the elliptical normal joint distribution 
with data, (Table [1] and Fig. H]), it is worse than a fit by the Weibull distribution. 
Despite the fact taking account of anisotropy of wind components modifies the tail 
of the distributions, the Gaussian approximation is too crude to model wind speed 
components. However, using Gaussian distributions is a powerful approach because 
analytic calculation can be performed. 

3.2. The Rice distribution 

As a second order approximation and for the sake of simplicity, we now consider the case 
of elliptic distributions with non zero mean p u and fi v and equal variances a' u = a' v = o. 
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equal variances for the two wind components. Thus, we obtain 

M ( M 2 + /j, 2 \ T fM/2\ 
Pm(M; fi,a ) = — exp I ——\i (-—\ j ( 5 ) 



with /i = \jy? u + (/%. The distribution in Eq. (JSJ) is called the Rice distribution and Iq is 
the th order modified Bessel function of the first kind. 



3.3. The Rayleigh-Rice distribution 

The key assumption we do in the following is that, in the Rhone valley area, there 
are two regimes: (i) random flow: the wind components have zero means and same 
variances, and the Rayleigh distribution describes well the wind speed statistics; (ii) 
channeled flow: the wind components have different means, and the Rice distribution 
describes well the wind speed statistics. The resulting distribution is a sum of the 
distribution (i) and (ii) conditioned to the absence and presence of the channeled flow 
occurrence. In other words, in the long-term, the wind system is described as the 
superposition of different local dynamics at different intervals that has been coined by 
Beck and Cohen (2003) as superstatistics or statistics of statistics. This yields, for one 
wind component, for instance the zonal wind, u (similar expression would arise for the 
meridional component, v): 

p(u) = J aMp(u)d». (6) 

In the case of the two regimes scenario, a u {^) is bimodal, 

a u (jj) = (1 - a) 5(ii) + a5(fi - fi u ), (7) 




Figure 4. Distributions of wind speed (blue) in m s _1 and a fit by a Weibull function 
(dashed black), Elliptic function (dash-dot red), and weighted Rayleigh-Rice function 
(dotted magenta) from observations made at (a) Montelimar and (b) Orange. 
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where 6 is the Dirac function and a is here the weight corresponding to channeled flow 
events occurrence. The PDF of M = \/u 2 + v 2 , is thus a combination of a Rayleigh 
distribution (corresponding to li u = li v = 0) and a Rice distribution. Such a method 
has been used to include in wind statistics modeling the contribution of zero wind speed 
(Tackle and Brown 1978, Tuller and Brett 1983). As a result we obtain the Rayleigh- Rice 
distribution: 

p M (M;/i,a 2 ,a) = ^exp f~^J ( a + i 1 ~ «) ex P (~^J J ° (~^)) ' (8) 

where Z = / pm{M; li, a 2 , a)dM is the normalization factor. 

When comparing to data, the weighted Rayleigh-Rice function provides a much 
better fit (Fig. [4]) and estimates much better the wind power (Table [1]) than the Weibull 
distribution. We indeed obtain less than 1.5% relative error for all locations and even 
less than 0.3% in Montelimar. Table [TJ shows that the weighted Rayleigh-Rice function 
gives by far the best estimate of the wind power at Lyon, Montelimar and Orange. The 
values of a correspond very accurately to the probability of occurrence of the dominating 
mistral situations (Orieux and Pouget 1984). Thus, taking account of the existence of 
regimes in regions where mountain environments is a fruitful approach to quantify wind 
speed statistics and to estimate wind power potential. 

4. Non-Gaussian distributions 

4-1. Super- statistics 

Figure |3] shows that at SIRTA, where the environment is much less complex than at 
the 3 other weather stations, the PDFs of the two wind components are not accurately 
predicted by a normal distribution and displays heavy tails. Such heavy tails illustrate 
the departure of the wind component PDFs from the normal distribution model. We 
show here that super-statistics can help for the understanding the non-Gaussian shape. 
Since both wind components at the SIRTA station have very close statistical properties 
(Fig. [3k), for the sake of concision, we will focus on one component, u, throughout this 
section. 

The large tails of the wind component distributions at SIRTA (Fig. [3k) originate 
in the transient nature of the wind. Meteorological conditions indeed change on several 
time-scales (day, anticyclones duration, season). For instance, the strongest wind being 
typically recorded in the winter, whereas in summer, long-lasting anticyclonic conditions 
give lower wind speeds. This induces a change of the statistical properties (e.g. wind 
component variance) on several entangled time scales. Here, we model this by assuming 
Gaussian shape of wind distribution, but with changing standard deviation, a. The 
super-statistics of the wind components can be derived as follows (Rizzo and Rapisarda 
2004): 

p(u) = J f(/3)p(u;/3)d/3, (9) 
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where / is the probability distribution of the fluctuating variable (3 = 1/cr 2 . Figure^ 
shows the distributions of short samples of wind component speed u measured at the 
SIRTA observatory, rescaled by their standard deviation, a, computed over 4 hr time 
periods. We see that the distributions of wind component velocity acquired during short 
time samples is well rescaled by its standard deviation. We can therefore assume the 




Figure 5. (a): Distributions of short samples of wind component speed u, rescaled 
by their standard deviation a computed over 4 hr time periods. Black line indicates 
a Gaussian function of standard deviation 1. (b): Distribution of each sample 
standard deviation f3 = a~ 2 and its fit by a Gamma distribution, (c): Supcrstatistical 
distribution corresponding to Eq. (|12|) . Color codes for varying parameters m and n. 
(d): Wind speed component distributions from the SIRTA observatory measurements 
(O) u, and x, v) and their fit by the superstatistical distribution corresponding to 
Eq. (fT!?]) with parameters, m = 9, n = 0.019 (black dashed line). 
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wind component distributions to be 

p(u; /3) = J - exp (-/3u 2 ) . (10) 



Figure [5b shows the distribution of each sample standard deviation = a~ 2 and its fit 
by a Gamma distribution. We see that the distribution of (3 is consistent with a Gamma 
distribution, of parameters m and n: 

I [m + ljn \n I y n J 

where T is the Gamma function. Hence, combining it to Eq. (jUJ), we obtain 



M w ) = TT ( 1 + ^ ) > C 12 ) 



where S is the Euler integral of the first kind. This distribution is called the generalized 
Boltzmann factor. It has been obtained when attempting to generalize the entropy 
definition (Tsallis 1988, Beck and Cohen 2003). The distributions modeled with Eq. (fT2"|) 
display heavy power law tails reminiscent of rare events (Fig. |5b). Unfortunately, our 
data has not enough statistics to observe such rare events and we can not compare 
the tails of the distribution (Fig. [5)1) • Yet, note that the central part of the observed 
distributions matches very well the distributions of Eq. ( Tl2|) (Fig. [5jl). 

We now turn to wind speed distribution, which is at stake for wind energy. For the 
sake of simplicity, we assume that the components u and v are independent statistically. 
Besides, we had previously observed that: (i) there is no correlation between u and v 
(Fig. Eh- inset); (ii) u and v have similar statistics, i.e. both components are described 
by the same m and n. Therefore the wind speed distribution is obtained by computing 
the joint distribution in radial coordinates and integrated over the angle. Therefore, we 
obtain 

nM A L , 1 n 2 M 2 . , , 

p(M) = -. wT-Bf-, - )F m, -, 1, , 13 

^ V ; 2B(±, m-1) 2 v 2'2 ; v '2' ' 1 + nM 4 ^' v ; 

where B is the Euler integral of the first kind and F is the hypergeometric function. 
Again, this yields a heavy tail distribution (Fig. EK), but for the set of parameters 
that fits the data the most (Fig. |6b), the wide tails occur for very low probability, and 
Eq. ( TlBl fits very well the data. 

4-2. Physical origin of super- statistics: stationary solutions of stochastic differential 
equations 

For a more in-depth analysis of the physical processes, several studies model the wind 
speed or its components using stochastic differential equation (SDE) (e.g. Nfaoui et al. 
2004, Bernardin et al. 2009). A first order SDE for stochastic process u can be written 
as: 

^ = C («) + £ ( 14 ) 
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Figure 6. (a): Superstatistical distribution corresponding to Eq. (|13|) . Color 
codes for varying m and n parameters, (b): Wind speed distribution from the 
SIRTA observatory measurements (O, w j and x, v) and its fit by the superstatistical 
distribution given by Eq. (|13|) with parameters m — 9, n — 0.019 (magenta dashed 
line) . 



where ( and £ are arbitrary functions and rj is a random function of time, often referred 
to as a "noise term". Considering the time evolution of the probability density function, 
p, associated to the stochastic variable, u, Eq. ( I14jl can be converted into the Fokker- 
Planck equation, 

£~S<« + 5h(<S«4 (15) 

If £ is constant, the system is said to be subject to additive noise, otherwise it is said 
to be subject to multiplicative noise. The stochastic processes u is usually a Markov 
process, i.e. characterized as memoryless: the next state depends only on the current 
state and not on the sequence of events that preceded it. A number of authors have 
modeled explicitly the temporal dependence of the wind speed through Markov chain 
models (Youcef Ettoumi et al. 2003, Nfaoui et al. 2004, Shamshad et al. 2005, Fawcett 
and Walshaw 2006a, Fawcett and Walshaw 2006b). Interestingly, it turns out that the 
distribution of Eq. ffl~2]) is also the stationary solution of Eq. f fT5|) . with C{u) = —'-/u and 
£(tt) = \/2jm\u\ + a/2x4, where and 7a are constants governing the strength of the 
multiplicative (first term) and additive (second term) noise. The relations between jm, 
7a and the parameters of Eq. (fTUI) are: 

7 + 7m 
m = 

7 7 ^ ( 16 ) 
n = 

1A 
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The multiplicative and additive noise terms can be physically interpreted as an interplay 
between turbulence, chaotic atmospheric variability and the mean wind (Bernardin 
et al. 2009). In the case of non-multiplicative noise, 7a/ = 0, the strength of the 
noise does not depend on the mean wind and the resulting stationary distribution is 
Gaussian. Conversely, the multiplicative stochastic process, used to model wind stress 
over the ocean by Sura and Gille (2003), makes the link between the wind variability 
and its mean state and induces a stationary power-law distribution (Sakaguchi 2001). 
Yet, as the observed wind component distributions are narrower than power-laws, 
this multiplicative stochastic model is crude to describe them. Finally, besides the 
rather accurate but empirical fit of the Weibull distribution to measured wind speed 
distribution, is there thus a physical reason to model wind speed statistics with a Weibull 
distribution? In order to address this question, a first step might be to consider a SDE 
with ((u) = — 7« and £(w) = y/2rfM\u\ k (Queiros 2008). Indeed, this yields the joint 
PDF for wind speed, M and direction, 0, 



with Z the corresponding normalization factor. Despite the fact that Eq. (fT7|) can not 
be integrated with respect to between and 2tt because of integration singularities, 
it shares common features with a Weibull distribution, namely the power-law and 
stretched-exponential with exponent k! = 2(1 — k). Future works will thus explore 
in more detail the potential of such SDE to model wind speed statistics. 

5. Conclusion 

The use of the Weibull distribution for wind statistics modeling and wind energy 
evaluation is a convenient and powerful approach. Is is however based on empirical 
more than physical justification and might display at times strong limitations for its 
application. In this article, we model the wind components instead of the wind speed 
itself, and thereby provide more physical insights on the validity domain of the Weibull 
distribution as a possible relevant model for wind statistics. 

Based on wind measurements collected nearby Paris in a rather flat environment 
and in Southern France in the Rhone valley in a ridge and mountains environment, this 
article proposes different probability density functions of the wind components. This 
opens the path to the analysis of the applicability of the Weibull distribution as a suited 
model for wind statistics and the quantification of the errors on the mean wind speed 
and wind power arising from the use of the Weibull distribution. 

This study first shows that in complex terrain, the existence of terrain induced wind 
system affects the surface wind distribution which does not display a Weibull shape. 
Modeling the wind statistics with a Weibull distribution produces errors exceeding 10% 
on the wind power estimate. The combination of a Rayleigh distribution to model the 




(17) 
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"isotropic" wind and a Rice distribution to model persistent wind regimes gives excellent 
results for the three surface weather stations located in the Rhone valley. Combining 
Rayleigh and Rice distributions is equivalent to applying the recently developed super- 
statistics theory, which models the wind system as the superposition of different local 
dynamics at different intervals with different mean wind speeds. However, Rayleigh 
and Rice distributions assume normal distributions for the two wind components. The 
distributions of the wind components, even over flat terrain, however exhibit non 
Gaussian distribution. This can also be modeled using super-statistics theory, which 
assumes fluctuating variances of the two wind components, that are eventually modeled 
by a Gaussian distribution over a short time interval. However, this produces power- 
law distribution for large wind speed value which contradict observed distributions and 
Weibull model. This can be partly overcome by assuming the wind components to verify 
stochastic differential equation with multiplicative noise with power coefficient. 

Parametric statistical downscaling using the Weibull distribution is often used to 
produce regional surface wind speed climatologies (Pryor et al. 2005). Also, this study 
points out the limits of the use of parametric distribution to downscale surface wind 
speed and advocates for non parametric statistical models (e.g. Michelangeli et al. 2009, 
Salameh et al. 2009, Lavaysse et al. 2012,Vrac et al. 2012). 
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